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ABSTRACT 



Context. Evershed clouds (ECs) represent the most conspicuous variation of the Evershed flow in sunspot penumbrae. 
Aims. We determine the physical properties of ECs from high spatial and temporal resolution spectropolarimetric measurements. 
Methods. The Stokes profiles of four visible and three infrared spectral lines are subject to inversions based on simple one-component 
models as well as more sophisticated realizations of penumbral flux tubes embedded in a static ambient field (uncombed models). 
Results. According to the one-component inversions, the EC phenomenon can be understood as a perturbation of the magnetic and 
dynamic configuration of the penumbral filaments along which these structures move. The uncombed inversions, on the other hand, 
suggest that ECs are the result of enhancements in the visibility of penumbral flux tubes. We conjecture that the enhancements are 
caused by a perturbation of the thermodynamic properties of the tubes, rather than by changes in the vector magnetic field. The 
feasibility of this mechanism is investigated performing numerical experiments of thick penumbral tubes in mechanical equilibrium 
with a background field. 

Conclusions. While the one-component inversions confirm many of the properties indicated by a simple line parameter analysis 
(Paper I of this series), we tend to give more credit to the results of the uncombed inversions because they take into account, at least 
in an approximate manner, the fine structure of the penumbra. 
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1. Introduction 

Evershed clouds (ECs) are patches of enhanced Doppler signals 
, that move in the penumbra of sunspots. Surprising ly, there ex- 

■ ist ver y few measurements of their magnetic fields. IShine et aP 
[ (|1994|) followed the evolution of ECs using Dopplergrams 
. and longitudinal magnetograms taken at the Swedish Vacuum 

■ Telescope on La Palma. They observed cases of ECs with weak 
] magnetogram signals, but this result could not be interpreted un- 
« ambiguously because the magnetograms were sensitive to both 
' field strength and field inclination variations. 

ICabrera Solana et aP (1200 7. hereafter Paper I) used simple 
parameters derived from visible and infrared lines to character- 
ize the polarimetric properties of ECs. While a complete obser- 
vational description was possible, the physical mechanism be- 
hind the EC phenomenon could not be studied without a more 
involved analysis of the measurements. 

Here we determine the physical properties of ECs by means 
of Stokes inversions, in an attempt to shed light on their na- 
ture. Two diflTerent models are employed: one-component atmo- 
spheres, and uncombed models that take into account the fine 
structure of the penumbra (Sect. [3]). Sections|4]and[5]describe the 
physical properties of the ECs and the penumbral filaments host- 
ing them, as well as their variation across the penumbra. Based 
on these results, in Sect. [6] we examine diflTerent mechanisms that 
could explain the EC phenomenon. 
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2. Observations 

Eor a complete description of the data sets the reader is referred 
to Paper L Time sequences of a small portion of the center- side 
penumbra of AR 10781 were taken with the two polarimeters of 
the German Vacuum Tower Telescope (VTT) on Tenerife. The 
polarimeters were operated simultaneously to record the Stokes 
parameters of four lines at 630 nm and three lines at 1565 nm. 
The cadence of the scans was 3.9 min on June 30, 2005 and 2 
min on July 1, 2005, when the spot was located at heliocentric 
angles of 43° and 35°, respectively. Eifteen ECs moving from 
the mid to the outer penumbra were identified as structures of 
enhanced Doppler velocities during the 236 min covered by the 
scans. Thanks to the adaptive optics system of the VTT, the spa- 
tial resolution of the infrared measurements (-Off 6) is one of the 
highest ever reached in ground-based spectropolarimetry. 



3. Stol^es inversions 

3.1. One-component models 

To invert the observed Stokes profiles we first use one- 
component models, i.e., we assume the whole resolution el- 
ement to be occupied by a single, laterally homogeneous at- 
mosphere. It is believed that penumbrae are formed by at 
least two magneti c components with different fi eld inclinations 
(Title et al." 1993; 'Solanki & Montavon' '1993"; "Martinez Pille^ 
2000; Bellot Rubio et al. 2004; Bello Gonzalez et al. 2005: 
Langhans et al.ll20Q5l: iBorrero et al .1120051: ISchlichenmaier et afl 
20071) . so the physical quantities retrieved from these inversions 
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Fig. 1. Illustration of the atmospheres filling the resolution ele- 
ment in an uncombed penumbra. The areas enclosed by the red 
and blue lines represent the background and tube atmospheres, 
respectively. The embedded tube is characterized by its center 
position (Ct) and radius (Rt). Ut and Lt are the upper and lower 
boundaries of the tube. The j-axis is the sunspot's radial direc- 
tion and z the local vertical. 



must be understood as an average of the various magnetic at- 
mosp heres coexisting in the resolution element (|Bellot Rubiol 
|2QQ3|) . The simplicity of the model may lead to erroneous con- 
clusions, but it is instructive to perform such inversions to allow 
comparisons with earlier works. 

The temperature stratification of the atmosphere is modified 
with 3 node 5j. As demonstrated in Paper I, visible and IR lines 
display diff'erent Doppler velocities. Hence, we allow for a gra- 
dient of LOS velocity (i.e., 2 nodes) to fit all the lines simulta- 
neously. In contrast, the magnetic field is supposed to be con- 
stant along the atmosphere. The inversion also returns height- 
independent macro- and micro-turbulent velocities, modeling 
the presence of unresolved motions. As the amount of stray light 
contamination is expected to depend on wavelength, the inver- 
sion code determines this contribution separately for each spec- 
tral range, which adds two more free parameters. In this way, the 
Stokes profiles emerging from the pixel are computed as: 



/ = (1 - a)Im + al, 



stray ? 



(1) 



where a is the stray light factor, /stray the average quiet Sun inten- 
sity profile representing the stray light contribution, and /m the 
Stokes vector emerging from t he model atmosphere . Our one- 
component models and those of iBorrero et al.l (|2004|) are diff'er- 
ent, since we consider height-independent magnetic fields. 

The inversions are carried out using the SIR code 

dRuiz Cobo & del Toro Iniestall 19921) with a total of 12 free pa- 
rameters. 



3.2. Uncombed models 

The typical size of the fine structure of the penumbra is 150- 
250 km or even smaller. At the spatial resolution of our obser- 
vations, atmospheres with very diff'erent properties are likely to 
be mixed in the resolution element. These atmospheres must 
be interlaced not only in the horizontal but also in the ver- 
tical direction to explain the non-zero Stokes V area asym- 
metries observed in visible and infrared lines. The uncombed 
model proposed by Solanki & Montavon (1993) incorporates 



^ Nodes are grid points of the model atmosphere where perturbations 
are sought. The new model atmosphere is constructed by interpolating 
the perturbations at the nodes to all grid points. 
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Fig. 2. Example of an initial uncombed model consisting of the 
background (dashed) and tube (solid) atmospheres. The gaus- 
sian is initially placed at logr = -0.8 and has a half-width of 
Alogr = 0.5. 



the interlacing of field lines in a natural way. This model 
envisages the penumbra as a collection of horizontal flux 
tubes embedded in a more vertical background field, which 
agrees with the results of spectropolarimetric analyses (e.g. , 
del Toro Iniesta et al. 2001; Schlichenmaier & Collados "2002!; 
Bellot Rubio et al. 2004; Borrero et al. 2004, 2005, 2006; Becl^ 
120061) . 

To account for the uncombed structure of the penumbra, at 
least to first order, the horizontal and vertical interlacing of the 
flux tubes and the background must be considered. To this end, 
we use a two-component model in which a fraction of the pixel is 
assumed to be filled by a background atmosphere and the rest is 
occupied by the background and tube atmospheres interlaced in 
the vertical direction (cf. Fig. [TJ. The emergent Stokes profiles 
are then the combination of the background, /b, and tube, /t, 
profiles: 



/=(l-a)[//t + (l-/)/b]+ a/stray. 



(2) 



The filling factor / represents the fractional area of the pixel 
occupied by the flux tube. 

The tube atmosphere is constructed adding gaussian pertur- 
bations to the stratifications of the background component (see 
Fig. [2]). The half-width (Rt) and position (Q) of the gaussian 
are free parameters derived from the inversion and are the same 
for all the atmospheric parameters. In contrast, the amplitudes 
of the gaussian depend on the parameter. For obvious reasons, 
we identify Rt and Q with the radius and central position of the 
tubes. Except for the temperature, the atmospheric quantities of 
the backgroun d are assumed to be constan t with optical depth 
The results of Idel Toro Iniesta et al.l (l200ll) . iBeflot Rubio et al.1 
(2004), and i Borrero et aP (l2004 demonstrate that the back- 
ground atmosphere harbors much smaller velocities than the 
tube component, so we take the background to be at rest. 

The uncom bed inversions have been performed with the 
SIRGAUS code (iBeUot Rubiol2003h . This code does not account 
for diff'erent stray light contaminations in the diff'erent spectral 
ranges, but systematic errors in the determination of the stray 
light should produce little changes of the magnetic field proper- 
ties. The number of free parameters is 17 (Table [TJ. 
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Table 1. Number of nodes employed in the uncombed inversion. 
- 1 indicates that the two atmospheric components have the same 
value of the corresponding parameter. 



Parameter 


Background 


Tube 


temperature 


2 


1 


velocity 





1 


field strength 


1 


1 


inclination 


1 


1 


azimuth 


1 


1 


microturbulence 


1 


1 


macroturbulence 


-1 


1 


stray light 


-1 


1 


filling factor 


-1 


1 


gaussian position 




1 


gaussian width 




1 






V 













3.3. Quality of the fits 

As an example, Fig. [3] shows the best-fit profiles resulting from 
the one-component inversion of a pixel located at a normalized 
radial distance r = 0.3 (r is defined such that represents the 
inner penumbral boundary and 1 the edge of the spot). The ob- 
served profiles are rather asymmetric, with Stokes U exhibit- 
ing four lobes. Taking into account the simplicity of the model, 
the match between observed and synthetic profiles is reason- 
ably good in both spectral ranges. This suggests that the scenario 
adopted for the one-component inversions is appropriate to de- 
termine the 'average' properties of the atmosphere. 

Figure m displays the same observed spectra and the best-fit 
profiles resulting from the uncombed inversion. The Stokes V 
area asymmetries and the shapes of Stokes Q and U are now re- 
produced in much greater detail. The use of the uncombed model 
considerably improves the quality of the fits as a result of its 
more realistic assumptions. 

4. One-component view of the EC phenomenon 

4.1. Physical properties of ECs 

Figures [5] and ?? show the atmospheric parameters inferred 
from the inversion of the June 30 data set. The diflTerent pan- 
els display the magnetic field strength, magnetic field inclina- 
tion in the local reference frame (LRF), the temperature T at 
logr = 0, the LOS velocity vlos at logr = -0.5, and the diflTer- 
ence between the unsigned LOS velocity in two different layers, 
Avlos = |vLos(logT = -0.5)1 - |vLos(logT = 0)|. This last quan- 
tity is proportional to the gradient of velocity with optical depth. 

The vector magnetic field displays the typical behavior ob- 
served in sunspots: it is weaker and more horizontal as the radial 
distance increases. It does not drop to zero outside the visible 
border of the spot, which is the signature of the sunspot mag- 
netic ca nopy feolank i et aLlll994l:lM^stendorp Plaza et al .1120011: 
iRezaei et al. 2006) . The ECs are visible as coherent structures 
with increased LOS velocity and field inclination. The pres- 
ence of many moving magnetic features in the surroundings of 
the spot is obvious (see the arrows plotted in the first panel of 
Fig. ??). Their origin and nature will be discussed elsewhere 
(ICabrera Solana et al ]|200i. 

For each EC, we calculate the maximum field strength, field 
inclination in the LRF, and LOS velocity at logr = -0.5 fol- 
lowing the procedure described in Sect. 6 of Paper L Assuming 
that the velocity and magnetic field vectors are aligned, the 
modulus of the velocity vector (the flow velocity) is simply 




Fig. 3. Example of Stokes profiles observed at r = 0.3 (dots) 
and best-fit profiles resulting from the one-component inversions 
(solid lines). The O2 teluric blends present in the visible spectra 
were not inverted, but they appear in the best-fit profiles due to 
stray light contamination. The vertical solid lines indicate the 
rest wavelengths of the lines. The Stokes V signal around Tii 
630.38 nm has been multiplied by six for better visibility (to the 
right of the vertical dotted line). 
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Fig. 4. Same as Fig.[3l for the uncombed inversion. 

M = Vlos/ cos 7los- We also determine the maximum differ- 
ence between the flow velocities at logr = -0.5 and logr = 
(A|v|). The results are summarized in Table [31 

The maximum field strengths attained by the ECs range from 
930 G to 1070 G with an average of 1005 G on June 30, and 
from 1300 G to 1500 G with an average of 1408 G on July 1. In 
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Fig. 5. Top to bottom: Magnetic field strength, field inclination, temperature at logr = 0, LOS velocity at logr = -0.5, and 
diff'erences between the unsigned LOS velocities at logr = and logr = -0.5. White contours outline the ECs and the inner and 
outer penumbral boundaries. The letters at the bottom of the last panel label each EC. 



Paper I we discussed the compatibility of the linear-to-circular 
polarization ratios and LOS velocities with an scenario where 
the magnetic and velocity fields of the ECs are nearly horizontal 
to the solar surface. Table |2] confirms the existence of very in- 
clined magnetic fields in the ECs: the average field inclination is 
deduced to be of the order of 82° for the June 30 observations 
and 71° for the July 1 data set. 

Both |v| and A|v| display higher values on June 30, i.e. when 
the spot is farther from disk center. The values of |v| are around 
4.5 km s"^ and 2.6 km s"^ on June 30 and July 1, while A|v| 
is about -0.8 km s"^and -0.3 km s"^on June 30 and July 1, 
respectively. The changes in |v| and A|v| cannot be explained by 
LOS eff'ects. Thus, the properties of the spot seem to have varied 
between the two observations. 



4.2. Radial variation of the physical properties of ECs 

The variation with radial distance of the atmospheric parameters 
derived from the one-component inversion is studied in this sec- 
tion. We calculate the radial curves of Fig. [6l in the same way 
as described in Paper I. The comparison of the physical proper- 



ties of the ECs and the penumbral filaments hosting them (intra- 
spines, according to Paper I) shows that: 



- On average, the magnetic field of the ECs is weaker than that 
of the intra- spines from the inner to the mid penumbra, but 
stronger in the outer penumbra. 

- ECs harbor more inclined fields than the intra-spines at all 
radial distances. 

- The azimuth of the field in the ECs is very similar to that 
found in the intra-spines at all radial distances, except per- 
haps in the outer penumbra. 

- The flow velocity inside the ECs is greater than in the intra- 
spines at all radial distances. 

- ECs display stronger gradients of the flow velocity with op- 
tical depth than the intra-spines. 

- ECs are slightly hotter than the intra-spines in the inner 
penumbra, but cooler from the mid to the outer penumbra. 



Therefore, in terms of one-component models, the magnetic, 
dynamic and thermal properties of the ECs are rather difl'erent 
from those of the intra-spines along which they move. 
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Fig. 6. Radial variation of the magnetic field strength, inclination in the LRF, azimuth, flow velocity, A|v|, and temperature for the 
ECs (solid lines) and the intra-spines along which they move (dashed lines), as inferred from the one-component inversions. Shaded 
areas indicate the rms fluctuations of the parameters, r = represents the inner penumbral boundary and r = 1 the outer edge of the 
penumbra 



Table 2. Physical properties of the observed ECs: field strength, 
inclination in the LRF, flow velocity at logr = -0.5, and the 
difl'erence between the flow velocities at logr = -0.5 and 
log T = 0. Errors represent the rms fluctuations of the physical 
parameters within each EC. 



EC 


B 


r 


|v| 


A|v| 




[G] 


[deg] 


[km/s] 


[km/s] 


A 


1020 + 20 


81+3 


4.5 + 0.4 


-1.0 + 0.6 


B 


1010 + 20 


83 + 2 


5.0 + 0.5 


-0.6 + 0.4 


C 


1050 + 60 


80 + 3 


4.5 + 0.3 


-0.8 + 0.6 


D 


1010 + 50 


82 + 5 


4.5 + 0.8 


-0.6 + 0.6 


E 


980 + 30 


81+3 


4.2 + 0.5 


-1.1 + 0.8 


F 


1070 + 60 


88 + 5 


4.8 + 0.8 


-0.4 + 0.4 


G 


1010 + 50 


84 + 2 


4.2 + 0.4 


-0.6 + 0.2 


H 


950 + 40 


83 + 1 


4.7 + 0.3 


-0.3 + 0.2 


I 


1000 + 50 


78 + 4 


4.1+0.3 


-1.0 + 0.3 


J 


930 + 40 


83 + 1 


4.7 + 0.3 


-0.9 + 0.3 


K 


1030 + 40 


79 + 2 


4.3 + 0.3 


-1.3 + 0.1 


Mean 


1005 


82 


4.5 


-0.8 


L 


1470 + 20 


60 + 3 


1.4 + 0.5 


-0.2 + 0.1 


M 


1440 + 70 


73 + 3 


2.5 + 0.5 


-0.2 + 0.1 


N 


1430 + 110 


77 + 2 


3.1+0.4 


-0.3 + 0.1 





1290 + 30 


75 + 4 


3.2 + 0.2 


-0.2 + 0.1 


Mean 


1408 


71 


2.6 


-0.2 



4.3. ECs as perturbations propagating along intra-spines 

In Paper I we found that ECs exhibit larger linear polarization 
signals than the filaments hosting them. This led us to suggest 



that they possess more inclined magnetic fields to the LOS. If the 
magnetic field and the velocity vectors are parallel, any increase 
in the LOS inclination would result in smaller LOS velocities. 
Thus, one has to conclude that the large Doppler signals asso- 
ciated with the ECs can only be produced by an increase in the 
modulus of the velocity vector, i.e., by stronger Evershed flows. 

The maps of field inclinations and velocities displayed in 
Fig. [7] for EC A demonstrate that this interpretation is correct. 
The strong linear polarization signals seen in the EC are in- 
deed produced by larger field inclinations to the LOS (upper 
left panel), which come from more horizontal fields (lower left 
panel). The higher LOS velocities detected in the EC are the ob- 
served signatures of stronger Evershed flows (lower right panel). 

During their journey to the outer penumbral boundary, the 
ECs modify the magnetic field of the intra-spines along which 
they move, making it more horizontal locally. After the passage 
of an EC, the intra- spine recovers its original field inclination. 
Figure [8] illustrates this behaviour for the case of EC A (the ar- 
rows indicate the orientation of the vector magnetic field). 

In view of these results, one could define the ECs as quasi- 
periodic perturbations of the magnetic and velocity fields of the 
penumbral filaments (intra-spines) along which they move. The 
main properties of these perturbations are the following: (a) they 
produce enhancements of y and |v|; (b) the amplitude of the per- 
turbation increases as it reaches larger radial distances; (c) when 
the EC starts to vanish, the amplitude of the perturbation is re- 
duced; and (d) the modification of y is perpendicular to the direc- 
tion of propagation of the perturbation (the vertical component 
of the EC proper motions is small compared with the horizontal 
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Fig. 7. Left: Maps of magnetic field inclination in the LOS reference frame {top) and the LRF (bottom). Right: Maps of LOS 
velocities at logr = -0.5 (top) and flow velocity at logr = -0.5 (bottom). White contours mark EC A, the intra-spine hosting 
it, and the boundaries of the penumbra. Black contours enclose pixels with y, 7los, and |v| greater than 80°, 40°, and 4.2 km s"^ 
respectively. Each tickmark in the _y-axis represents ^^ 




t = 3.9 



t = 7.8 t=11.7 t=15.6 




t = 19.5 



Fig. 8. Variation of y along the cuts shown in Fig. [7] as EC A propagates from the inner = 3.9 min) to the outer (t = 19.5 min) 
penumbra. The color scale ranges from 55° (black) to 85° (red). Each arrow has the same length and indicates the orientation of the 
magnetic field vector retrieved from the inversion. White contours outline EC A and the penumbral boundaries. 



component, cf. Paper I). Except for (c), the behavior described 
in this section is also observed in other type I ECs. 



5. Uncombed view of the EC phenomenon 

The spine/intra- spine organization of the penumbra deduced 
from one-component inversions can be interpreted as the az- 
imuthal variation of the filling factors of two components whose 
properties remain relativ ely constant at a given radial distance 
dBellot Rubio et al.ll2004l) . Thus, the question naturally arises as 
to whether the more inclined fields and larger flow velocities of 
the ECs indicated by the one-component inversions are real or 
an artifact of simplistic modeling. 

Figures [9] and ?? show the physical properties of the tube 
component derived from the uncombed inversions, for the June 
30 data set. They already make it clear that the properties of 
the tubes in the ECs and the intra-spines are similar, except for 
the filling factor. In other words: when the fine structure of the 
penumbra is considered, the ECs and the intra-spines do not 
seem to have diff'erent physical properties. This is in contrast 
with the results of one-component inversions (Figs.[5]and ??). 



5.1. Properties of the tube and background atmospheres 
inside the ECs 

Figure [TOl displays the radial variation of the field strength, LRF 
inclination, azimuth, flow velocity, temperature at log r = 0, and 
filling factor of the tube and background atmospheres in the ECs. 
Similar curves are ob tained from classical tw o-component inver- 
sions of the spectra ("Cabre ra Solanal '2007). the only diff'erence 
being the weaker background fields indicated by the uncombed 
inversion^!. 

As expected, the field is more inclined in the tubes 
than in the background at all r adial distances. In agr e ement 
with previous investigations ( Westendorp Plaza et al. 1997', 
2001; Schlichenm aier & Schmidt 20 00; Mathew et al. 2003i 
iBellot Rubio et al.ll2004 iBorrero et al.ll2004 12005L [2006h . field 
lines diving back to the solar interior (y^ > 90°) are found in the 
outer penumbra at the position of the ECs. The tube and back- 
ground fields show diff'erences in azimuth, but they are usually 
smaller than 20° . 

Figure [TOl demonstrates that the flow velocity increases 
monotonically with radial distance within the tubes, varying by 
some 0.7 km s"^ from the inner to the outer penumbra. The tube 
temperatures at logr = are smaller than those of the back- 



^ These weak fields may not be real but a a consequence of initializ- 
ing the inversion with low temperatures in the tubes. 
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Fig. 9. Top to bottom: Magnetic field strength, field inclination, temperature at log r = 0, LOS velocity, and filling factor of the tube 
atmosphere retrieved from the uncombed inversion. The parameters of the tubes are constructed taking the values at the tube center. 
White contours outline the ECs and the inner and outer penumbral boundaries. The letters at the bottom of the last panel label each 
EC. 



gr ound at all radial distances. T he same behavior has been found 
by lBorrero et al.l (l2QQ5l l2QQ6l) , but only in the outer penumbra. 
The tube filling factor increas es monotonica lly with radial dis- 
tance, in agreement with .Borrero et al. l (l2QQ5h . 

5.2. Position and width of tlie tubes 

Figure fTTI examines the center position, radius, and boundaries 
of the tubes deduced from the inversion. At all radial distances, 
the tube axes are located above log r ~ 0, with similar values 
in the ECs and the intra-spines. Within the scatter we do not de- 
tect diff'erences between the ECs and the intra-spines in terms 
of width or position of the tube's upper boundary. The lower 
boundary is always below the line forming region, so the tubes 
are necessarily optically thick and/or low lying. This implies that 
the spectral lines only sample their upper halfs. A similar result 
has been reported by Borrero et al.^ (.2006.) from uncombed in- 
versions of penumbral profiles. 

5.3. ECs as structures of increased filling factor and flow 
velocity 

Figure [12] shows maps of the filling factor of the tube (/) as de- 
rived from the uncombed inversions (only the case of EC A is 
considered; the other ECs behave in a similar way). The maps 
clearly demonstrate that ECs are regions of increased filling fac- 
tors. The same conclusion can be drawn from the left panel of 



Fig. [131 where we plot the radial variation of / for EC A and the 
rest of the intra- spine. The flow velocity is also higher in the EC 
(right panel of Fig.[T3]). 

The mean difl'erences between the filling factors (A/) and 
flow velocities (A|v^|) of the ECs and the intra-spines are dis- 
played in Fig. [14] as a function of radial distance. A/ is of the 
order of 10%, while A|v^| reaches 0.7 km s"^ Similar plots for 
the field strength, inclination, azimuth, and temperature diff'er- 
ences are also presented in Fig. [14] The values of and ij/i turn 
out to be remarkably similar in the ECs and the intra-spines. The 
same happens with and T^. This implies that the uncombed 
inversion no longer interprets the ECs as perturbations of the 
magnetic configuration of the intra-spines. Figure [15] shows the 
corresponding quantities for the background atmosphere: while 
the passage of ECs is not associated with changes of yb or j/^b, 
we observe small decreases of 5b and slight increases of Tb in 
the inner penumbra. 

To understand why the code returns enhanced filling factors 
in the ECs, Fig. [16] compares the Stokes profiles emerging from 
the same pixel before and during the passage of EC A (cf. the 
numbers marked in Fig. [12]). The figure also displays the contri- 
bution of the background and tube atmospheres to the profiles. 
Note that (a) the EC passage increases the linear-to-circular po- 
larization ratios mainly through an enhancement of the linear 
polarization signal, and (b) the linear polarization arises almost 
exclusively from the tube component. This combination of fac- 
tors allows the code to explain the EC profiles just by increasing 
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Fig. 10. Radial variation of the magnetic field strength, field in- 
clination, field azimuth, flow velocity, temperature at log r = 0, 
and filling factor of the tube (solid) and background (dashed) 
components within the ECs, as derived from the uncombed in- 
versions. Shaded areas represent the rms fluctuations of the at- 
mospheric quantities at a given distance. 
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Fig. 11. Center position, radius, and boundaries of the flux tubes 
in the ECs (solid lines) and the intra-spines hosting them (dashed 
lines). The quantities are given in units of the logarithm of the 
optical depth. The shaded areas represent the standard deviations 
of the parameters. 
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Fig. 13. Radial variation of the filling factor (left) and flow ve- 
locity {right) inside EC A (solid lines) and its intra-spine (dashed 
lines). 
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Fig. 14. Top: Difl'erences between the filling factor (left) and 
flow velocity (right) of the tubes inside the ECs and the intra- 
spines as derived from the uncombed inversions. The shaded 
areas represent the rms fluctuations of the difl'erences. Middle: 
Same, for the field strength and inclination. Bottom: Same, for 
field azimuth and temperature at log r = 0. 



Fig. 12. ]V[aps of the filling factor of the tube (/) derived from the 
uncombed inversions. White contours mark the intra-spine and 
the boundaries of the penumbra. EC A is outlined with blue con- 
tours. Black contours delimit pixels having / larger than 77%. 
Asterisks and numbers indicate the position of a pixel before 
and during the passage of the EC. 



the contribution of the tube atmosphere, i.e., /. In other words: 
the filling factors of the ECs are large as a consequence of the 
high linear- to-circular polarization ratios they exhibit (Paper I). 
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Fig. 15. Same as Fig. [14] for the field strength, inclination, az- 
imuth, and temperatures at log r = of the background. 
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Fig. 16. Observed linear and circular polarization profiles of Fe I 
1564.85 nm (dots) emerging from a pixel before (#1) and dur- 
ing (#2) the passage of EC A. The pixel is marked with asterisks 
and numbers in Fig. [121 Blue and green lines represent the con- 
tributions of the tube and background components to the best-fit 
profiles (red lines). 



5.4. Physical origin of the Increased filling factors 

The large filling factors of the ECs may result from (a) an in- 
crease in the number of penumbral tubes, and/or (b) a greater 
fraction of the tube's cross section in the line forming region. 

Option (a) seems unlikely because it would require the par- 
ticipation of new tubes, all of them migrating outward along 
intra-spines in an organized way as the ECs. Option (b) would 
require flux tubes with larger radii or placed higher in the photo- 
sphere. Both possibilities are ruled out by Fig.[TT] 



Therefore, we deem that perturbations of the thermodynamic 
properties of the tubes are more feasible as the source of the 
filling factor enhancements associated with the ECs and, con- 
sequently, as the origin of the EC phenomenon. Since no sig- 
nificant temperature variations occur, we suggest that the EC 
phenomenon is produced by variations of density/pressure in- 
side penumbral flux tubes. This mechanism will be examined in 
Sect.f 



6. Discussion 

6.1. One-component Interpretation 

In Paper I we found that ECs are penumbral structures charac- 
terized by larger Doppler shifts, stronger linear polarization sig- 
nals, and larger Stokes V area asymmetries than both spines and 
intra-spines. We suggested that these were observational signa- 
tures of stronger Evershed flows and more horizontal magnetic 
fields in the ECs as compared with the rest of the penumbra. Our 
one-component inversions seem to confirm this picture. 

Two diff'erent mechanisms have been proposed to explain the 
origin of ECs: magnet oacoustic waves superimposed on a s teady 
or quasi- steady flow (Georgakilas & Christopoulou '2003"), and 
Q or U kinks propagating al ong field lines (Ryutova et al. 1998i: 
iRouppe van der Voonl2003h . Another possibility is that ECs are 
the signatures of sea-serpent field lines associated with mov- 
ing penumbral tubes ( Schlichenmaier 2002). None of the three 
mechanisms seems to be compatible with the magnetic field ge- 
ometry deduced from the inversions. The reason is that they 
should produce both an increase and a decrease of the field in- 
clination in the ECs, but the former is not detected. 

We cannot exclude that higher resolution observations may 
solve the problem in the future. With the present data, however, 
the conclusion is that origin of the magnetic field perturbations 
indicated by the one-component inversions remains unknown. 

6.2. Uncombed Interpretation 

We have carried out simple numerical experiments to examine 
whether or not larger linear- to-circular polarization ratios and, 
consequently, an enhancement of /, may result from density and 
pressure variations inside penumbral flux tubes. To that end we 
use the model of flux tubes in mechanical equilibrium proposed 
by IBorrerd (12007). In this model the tubes are not thin, i.e., the 
physical properties change over their cross sections. 

We assume two difl'erent equilibrium configurations to 
model the conditions of the tube and background atmospheres in 
intra-spines before and during the passage of an EC. In order to 
emulate the conditions before the EC passage we consider a hor- 
izontal (rt = 90°) tube with a field strength of 1000 G that chan- 
nels an Evershed flow of |vt| = 4.5 km s"^ . The tube is located at 
the heliocentric angle of the observations (43° on 30 June) with 
its axis pointing along the line of symmetry. It is surrounded by 
a background field with yb = 60° and = 450 G. The upper 
panels of Fig. [T7] show the temperatures, gas pressures, and den- 
sities derived from the model. The x-axis is perpendicular to the 
tube axis and the z-axis represents the line of sight. 

The only property that seems to change during the passage 
of an EC is the strength of the background field. Thus, we model 
the new situation using the same parameters for the tube com- 
ponent (rt = 90°, Bt = 1000 G, rb = 60°) but a weaker field 
of B\y = 350 G in the background atmosphere. To maintain lat- 
eral pressure balance, the gas pressure and density of the tube 
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Fig. 17. Top panels: Temperature, gas pressure, and density in the tube and their surroundings simulating the conditions of the 
atmosphere before the EC passage. Bottom panels: Variation of the temperature, gas pressure, and density produced by the EC. 
The tube is enclosed by the black lines. Solid and dashed white lines indicate isocontours of optical depth (logr = -1, -0.5, 0, 0.5) 
before and during the EC passage, respectively. 



decrease in response to the weakening of the background field 
(bottom panels of Fig. [TT]) . 



The pressure and density variations associated with the EC 
move the optical depth scale toward deeper layers, which in- 
creases the fraction of the tube inside the line formation region 
(compare the solid and dashed lines in the bottom panels of 
Fig.[T7]). It is important to remark that the variation of the optical 
depth scale is negligible in the upper half of the tube, becoming 
significant only in the lower half. This may explain why the po- 
sition of the tube's upper boundary is not seen to change with 
the EC passage. 



The downward shift of the optical depth scale increases the 
observed L/V ratios. This is demonstrated in Fig. [181 where we 
show the total linear and circular polarization profiles of Fei 
1565 nm before and during the EC passage. The contributions 
of the tube and background atmospheres to the emergent pro- 
files are indicated with blue and green lines, respectively. When 
the EC moves along the intra-spine, the linear-to-circular polar- 
ization ratio grows from LjV = 0.22toL/y = 0.27. Since the 
amount of circular polarization remains roughly the same, the 
increase in L/V is due to an increase in the linear polarization, 
much in the same way as observed (Fig. [161). The larger LjV 
ratio is mainly due to a larger flux tube contribution. Since the 
magnetic properties of the tube are the same before and during 
the EC passage, the only possible cause is a change in the ther- 
modynamic properties of the tube: a decrease in the gas pressure 



and densit^ The larger flux tube contribution is interpreted by 
uncombed inversions as enhancements of the filling factor. 

Our experiments support the idea that the L/V increase ob- 
served during the passage of ECs is the result of lower gas 
pressures and densities in the penumbral tubes that carry the 
Evershed flow. With the model of Borrer^ (l2007l) , these condi- 
tions are forced imposing a weaker external magnetic field. Such 
a field is consistent with the observations (Fig. [151), but its origin 
remains unclear. 

7. Summary 

In this paper we have interpreted the Stokes profiles emerging 
from Evershed clouds (ECs) using simple one-component and 
more sophisticated uncombed inversions of 4 visible and 3 in- 
frared lines. The inversions have allowed us to determine the 
physical properties of ECs, the ultimate goal being to understand 
the nature of the EC phenomenon. 

The one-component inversions confirm many of the prop- 
erties indicated by the line parameter study of Paper I. The re- 
sults of these inversions suggest that ECs are structures having 
more inclined fields and stronger Evershed flows than the rest 
of the penumbra. ECs move along intra-spines toward the outer 
boundary of the spot. In their journey, they change the magnetic 
configuration of the intra-spines, which only recover their initial 

^ To confirm this point, we have repeated the experiment assuming 
that the gas pressure in the tube increases during the passage of the EC. 
This is induced through a larger background magnetic field (5b = 550 G 
vs 5b = 450 G). As expected, the larger gas pressures and densities in 
the tube shift the optical depth scale toward higher layers, reducing the 
contribution of the tube to the linear polarization profiles and, thus, the 
L/y ratio. 
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Fig. 18. Synthetic linear and circular polarization profiles of Fe i 
1564.8 nm (black curves) simulating the conditions before (top) 
and during (bottom) the EC passage. Blue and green profiles rep- 
resent the contributions of the tube and background atmospheres 
to the emergent profiles, respectively. 



state after the passage of the ECs. Thus we suggest that the EC 
phenomenon is caused by a perturbation of the magnetic field of 
the intra-spines. None of the scenarios proposed so far for the 
EC phenomenon seem to be capable of explaining these find- 
ings. Both magnetoacoustic waves and magnetic kinks should 
give rise to vector fields pointing to the solar interior, but we do 
not detect them. 

When we account for the fine structure of the penumbra, the 
EC phenomenon is no longer interpreted as a real perturbation of 
the magnetic field of the intra-spines. The uncombed inversions 
suggest that the physical properties of the penumbral flux tubes 
are rather similar in the ECs and the intra-spines hosting them. 
The background atmospheres are also similar, except for a small 
decrease in the field strength and slightly larger temperatures. 
The only significant diff'erence observed in the intra-spines dur- 
ing the passage of an EC is an enhancement of the filling factor 
of the tube component, i.e., a larger visibility of the flux tubes 
that carry the Evershed flow. The inversions indicate that the en- 
hanced visibility is not due to changes in the position or radius 
of the tubes. We therefore propose that the EC phenomenon is 
the result of variations of gas pressure and density in the flux 
tubes. Such variations would propagate toward the outer sunspot 
boundary, producing the motion of the ECs. Our interpretation 
is summarized in Fig. \19\ Type I ECs vanish at the edge of the 
penumbra because the tubes return to the solar interior and the 
density/pressure perturbations go out of the line forming region. 
Type II ECs disappear in the sunspot moat (cf. Paper I), so the 
flux tubes associated with th em have to continue well bey ond the 
outer penumbral border (see lCabrera Solana et al ] l20Q8l) . 

We have performed simple numerical experiments to exam- 
ine the feasibility of this scenario. We find that a decrease of 
the gas pressure in the upper half of the tube implies a reduc- 
tion of the density, which shifts the optical depth scale toward 
deeper atmospheric layers. The shift does not change the optical 
depth of the tube's upper boundary, but makes the tube occupy 
a larger fraction of the line forming region. The result is an in- 
crease in the linear-to-circular polarization ratios associated with 
the ECs, which are interpreted by the inversion code as larger 




Inner penumbral boundary 



Outer penumbral boundary 



Fig. 19. Cartoon of a type I EC propagating along a penum- 
bral flux tube which dives back to the solar interior at the outer 
sunspot edge. The EC is represented by the black circles. The 
background field is indicated by the arrows. The dashed line 
marks the r = 1 level, and the vertical dotted lines the penumbral 
boundaries. For visualization purposes, the spot is significantly 
compressed in the horizontal direction. This makes the flux tube 
appear to be elevated, but in reality it is a deep-lying tube that 
never reaches high photospheric layers. 



filling factors. This mechanism holds promise to explain the EC 
phenomenon. However, further numerical and theoretical work 
is needed to understand the origin of the pressure/density pertur- 
bations. An interesting possibility is that they are the signatures 
of shocks inside pen umbral flux tubes, like the ones detected by 
iBorrero etaD(l2005h . 
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